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Front cover: Three southern flying squirrels (Glaucomys volans) emerging from a nest cavity in a mature dead pine 
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Abstract 

Ecological surveys were conducted in central Mexico to assess parameters associated with 
the distribution, habitat preference, and nesting sites of the southern flying squirrel (Glaucomys 
volans). These surveys indicated that deforestation and harvesting of mature oak trees has de¬ 
populated flying squirrels from a large portion of their historic range. DNA sequences obtained 
from the mitochondrial control region were used to evaluate genetic divergence among four 
populations in central and southern Mexico. Levels of intrapopulation genetic variation were 
greatest in the Chapa de Mota, Mexico population (average genetic distance = 2.8%), whereas 
interpopulational comparisons showed the Chapa de Mota, Mexico population to be the most 
different from the Floresta, Chiapas population (average genetic distance = 4.1%). Results 
indicate that substantial genetic variation exists between populations in southern Mexico and 
suggests either a historical bottleneck associated with the Isthmus of Tehuantepec in Chiapas or 
more widespread gene flow across central Mexico influenced by the Sierra Madre Oriental and 
Trans-Mexican Volcanic Zone. Conservation efforts are needed, especially in southern Mexico, 
to preserve populations of southern flying squirrels threatened by habitat loss. 
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Resumen 

Se condujeron inventarios ecologicos en la parte central de Mexico, con el fin de evaluar 
los parametros asociados a la distribucion, preferencia de habitat, y sitos de crianza de la ardilla 
voladora surena (Glaucomys volans). Estos inventarios indicaron que la deforestacion y la 
recoleccion de encinos maduros han despoblado a las ardillas voladoras de una gran porcion 
de su distribucion historica. Secuencias de ADN obtenidas de la region control mitocondrial 
fueron utilizadas para evaluar la divergencia genetica entre cuatro poblaciones del centra y 
sur de Mexico. Los niveles de variacion genetica fueron mayores en Chapa de Mota, Mexico 
(distancia genetica media = 2.8%), mientras que las comparaciones interpoblacionales muestran 
que Chapa de Mota, Mexico es muy diferente que la poblacion de Floresta, Chiapas (distancia 
genetica media = 4.1%). Los resultados indican que existe una variacion genetica substancial 
entre poblaciones en el sur de Mexico y sugiere que esto podria ocurrir por cuellos de botella 
historicos asociados con el istmo de Tehuantepec en Chiapas o por flujo genico a lo largo del 
centra de Mexico influenciado por la Sierra Madre Oriental y la Zona Volcanica Transversal de 
Mexico. Esfuerzos de conservacion son necesarios, especialmente en el sur de Mexico, para 
preservar las poblaciones de ardilla voladora del sur amenazadas por perdida de habitat. 

Palabras clave: ardillas voladoras, conservacion, diversidad genetica, Glaucomys volans , 
Mexico 
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Introduction 


The southern flying squirrel ( Glaucomys volans ) 
is distributed from southeastern Canada through the 
eastern United States to central Honduras (Hall 1981). 
Typically, G. volans is associated with deciduous-hard- 
wood forests, relying on mature oak trees ( Quercus) for 
food and nest cavities (Dolan and Carter 1977). This 
reliance on mesophytic forests coupled with Pleistocene 
climatic changes has resulted in a disjunct distribution 
for G. volans (Braun 1988; Arbogast 1999; Ceballos 
et al. 2010). Three allopatric distributions have been 
documented, with the largest occurring in the eastern 
United States and southern Canada, the second occupy¬ 
ing southern Mexico and northern Central America, and 
the third including a small area in southwestern portion 
of the Mexican State of Chihuahua (Goldman 1936). 
Information on the present distribution and habitat use 
of flying squirrels in Mexico is sparse (Ceballos and 
Oliva 2005). The species is considered threatened in 
Mexico due to a restricted distribution and habitat loss 
(SEMARNAT 2002). 

Flying squirrels in Mexico and Central America 
probably originated from a population in the south¬ 
eastern United States, following corridors of suitable 
habitat during the Pleistocene, but the timing and 
dynamics of the interchange are unknown (Ceballos 
et al. 2010; Kerhoulas and Arbogast 2010). Arbogast 
(1999) suggested that the patterns of genetic varia¬ 
tion in the southern flying squirrel were the result of 
Pleistocene shifts in the distribution of oak ( Quercus ) 
and hickory ( Carya ) forests. Historically, six subspe¬ 
cies of flying squirrels were documented in Mexico 
and Central America ( chontali, goldmani, herreranus, 
madrensis, oaxacensis, and underwoodi. ; Hall 1981); 
however, only three subspecies currently are rec¬ 
ognized {goldmani, guerreroensis, and oaxacensis'. 


Diersingl980). Recently, Kerhoulas and Arbogast 
(2010) examined five populations from Mexico and 
Guatemala and reported that genetic divergence values 
in the mitochondrial cytochrome-# gene ranged from 
0.71% (Oaxaca Highlands and Sierra Madre de Chia¬ 
pas) to 4.38% (Chiapas Highlands and Sierra Madre 
del Sur). These five populations differed genetically 
from populations located in eastern North America by 
divergence values ranging from 2.3-4.9%. Kerhoulas 
and Arbogast (2010) concluded that two primary ge¬ 
netic groups existed - one from the Sierra Madre del 
Sur and another that included all other populations 
from Mexico and Guatemala, although the later group 
potentially could be subdivided into additional groups. 
These two groups were hypothesized to have diverged 
from one another approximately 0.75-0.50 million 
years ago. No definitive taxonomic conclusions could 
be reached because of incomplete sampling of some 
subspecies and geographic regions. 

In this paper, we evaluate the following issues 
regarding flying squirrels in Mexico: 1) geographic 
distribution and habitat use; 2) genetic diversity; and 
3) conservation status. We specifically assessed the 
genetic diversity of four populations of G. volans from 
the states of Chiapas, Mexico, and Queretaro. These 
populations represent different geographic regions, 
including the Sierra Madre Oriental and the Trans- 
Mexican Belt, which were not available for study by 
Kerhoulas and Arbogast (2010). Nucleotide sequences 
from the mitochondrial DNA control region (D-loop) 
were used to assess genetic divergence within and 
among populations. The D-loop region evolves at a 
sufficient rate (Petersen and Stewart 2006) to detect 
genetic divergence among populations in flying squir¬ 
rels. 


Methods 


Geographic distribution and habitat use. —The 
historical distribution of G. volans in Mexico was 
determined on the basis of published and unpublished 
records. Localities where species have been recorded 
were characterized for geographic coordinates, veg¬ 
etation type, and elevation. Localities were grouped 


together in eight more or less discrete mountain regions. 
The area of each region and the distance between clos¬ 
est regions were calculated. 

To determine habitat requirements of flying 
squirrels, we compared four adjacent areas in Pinal de 
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Amoles, Queretaro. This locality possesses a rough 
topography, covered by pine, oak and pine-oak forests 
surrounded by an arid zone. Vegetation associated with 
the zone has been highly fragmented due to agricultural 
and forestry practices. We randomly selected 12 plots 
(100 x 50 m) known to harbor flying squirrels and four 
plots without squirrels. In each plot, tree diversity, den¬ 
sity, height, and DBH (diameter at breast height) of all 
the trees > 10 cm of DBH were measured. Each time a 
nest was located we recorded the following variables: 
tree species, DBH, tree height, number of dens per 
tree, diameter and height of den, height and width of 
the entrance of the den, content of the den (e.g., feath¬ 
ers, moss, hair, lichen), and number of animals found 
inside the den. Tree diversity was calculated using the 
Shannon-Wiener diversity index (Zar 1984). 


Stepwise discriminant functions analysis (DFA) 
was used to separate plots that were occupied by 
squirrels from those without squirrels and to identify 
variables contributing to the distribution of squirrels 
in the zone. Prior to the multivariate analysis and 
univariate screen procedure (ANOVA), variables were 
selected having an F value with at least p < 0.10. An 
F-to-enter of 3.84 and an F-to-remove of 2.71 were 
utilized. Distribution of squirrels was investigated us¬ 
ing a DFA analysis, habitat variable, plots possessing 
squirrels (e.g., nests, squirrels present), and plots with 
no evidence of squirrels. 

Sampling .—Flying squirrels were collected from 
four localities in central and southern Mexico (Fig. 1) 
by either wiring a Sherman live-trap (H. B. Sherman 


102° 90° 



Figure 1. Geographic distribution of Glaucomys volans in Mexico (see Table 1 and Appendix I for details). Focalities 
(1-37) utilized in the habitat and nest cavity analyses are represented by boxes; in addition, localities 6, 12, 13, and 36 
were used in the genetic analyses and are indicated by circles. Roman numerals (I-VIII) reflect the major mountain 
ranges discussed in this study. 
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Trap Company, Tallahassee, Florida) to a tree limb 
(baited with bacon and apple bits, oats, or chicken 
meat) or by placing a net over a tree cavity and striking 
the tree and capturing individual squirrels in a hand 
net. Flying squirrels were anesthetized with ketamine 
(following methods approved in the guidelines of the 
American Society of Mammalogists Animal Care and 
Use Committee 1998) and hair, blood, or toe samples 
were collected and the animal released. In some cases, 
animals were euthanized and tissue samples (muscle, 
liver, heart, and kidney) were obtained. 

Genetic diversity .—Mitochondrial DNA was 
extracted either from hair, blood, or tissue samples 
and purified using the Wizard Miniprep kit (Promega, 
Madison, Wisconsin), or genomic DNA was isolated 
following the methods of Smith and Patton (1999). 
Depending on the quality of DNA, the entire D-loop 
region (approximately 1,000 bp) and a portion of the 
5’-proline tRNA(approximately 85 bp) were amplified 
using PCR primers specific to the D-loop (Castro- 
Campillo et al. 1999) or smaller regions of the D-loop 
(200--300 bp) were amplified using internal D-loop 
primers (Mendez-Harclerode et ah 2005). 

The following thermal profile was used for PCR 
amplification: an initial cycle of 93.5° C for 1 min, 
followed by 33 cycles of 93.5° C for 40 sec, 49° C for 
40 sec, and 72° C for 1 min 30 sec, and a final cycle 
of 72° C for 2 min. PCR products were purified with 
a QIAquick kit (Qiagen Inc., Chatsworth, California). 
Sequencing reactions used primers reported in Mendez- 
Harclerode et al. (2005). Cycle sequencing was con¬ 
ducted using the ABI Prism dRhodamine terminator 
ready reaction mix (Applied Biosystems, Foster City, 
California) or ABI Big Dye version 3.0 ready reaction 
mix and samples were analyzed on an ABI Prism 310 
or ABI 3100 -Avant automated sequencer (Applied 
Biosystems, Foster City, California). Both forward 
and reverse sequences were obtained. Sequences 
initially were aligned using Sequencher 3.1 software 
(Gene Codes, Ann Arbor, Michigan) and then adjusted 
manually. All sequences were deposited into GenBank 
(accession numbers are provided in Appendix I). 


Two samples of G. volans obtained from Gen- 
Bank (Petersen and Stewart 2006) were used as out¬ 
groups in all analyses. Variable nucleotide positions 
were treated as unordered, discrete characters with five 
possible states: A, C, G, T, or gaps. Several alignments, 
reflecting differences in the placement of gaps, were 
examined using the heuristic search option in PAUP* 
(Swofford 2002). The alignment that produced the 
shortest tree was assumed to have incorporated the best 
placement of gaps and was used in all analyses. 

The Tamura-Nei model of evolution (Tamura and 
Nei 1993), developed specifically for analysis of DNA 
regions with different rates of evolution, was used to 
calculate genetic distances. These values were then 
used to assess levels of genetic divergence within and 
among populations. A network was constructed using 
the software program TCS (version 1.18; Clement et 
al. 2000) under the method of statistical parsimony 
(Templeton et al. 1992). Twenty-two individuals 
(outgroup excluded) were included in this analysis and 
potentially unsampled haplotypes were estimated. 

Conservation Status .—To evaluate the pres¬ 
ent conservation status of flying squirrels in central 
Mexico, surveys were carried out in five localities in 
three states (Hidalgo, Queretaro, and Mexico). Lo¬ 
calities where relatively healthy populations of flying 
squirrels were found in the last decade were selected 
for this study (Ceballos and Galindo 1983; Ceballos 
and Miranda 1985). We searched for flying squirrels 
in pine-oak and oak forests using a simple searching 
technique that has been very useful to locate and capture 
flying squirrels (Ceballos and Galindo 1983; Ceballos 
and Miranda 1985). The technique consisted of estab¬ 
lishing transects in dense mature forests and searching 
all trees for crevices and nests. Once a hollow or nest 
was found, squirrels were captured with either a hand 
net or trapped as described previously. At each locality 
we recorded the number of nests and flying squirrels 
found, general conditions of the habitat, and the pres¬ 
ence of human activities impacting the forests. 
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Results 


Geographic distribution and habitat use. — Glau¬ 
comys volans has a scattered distribution in eastern, 
central, and southern Mexico (Fig. 1). Knowledge of 
the species’ distribution in Mexico has become more 
detailed in the last two decades, where most locali¬ 
ties along the Trans-Mexican Volcanic Belt in central 
Mexico have been recorded, including new records 
from Villa Victoria, San Cayetano, Chapa de Mota, and 
Dexcani (Mexico), Zitacuaro (Michoacan), and Santa 
Elena (Hidalgo). In Mexico, flying squirrels are known 
from 37 localities in 11 states (Fig. 1 and Table 1), with 
most records from the southern states of Oaxaca and 
Chiapas. The northernmost record is from Rancho 
del Cielo (Tamaulipas; 23°04' N, 99°IF W), southern¬ 
most from San Sebastian Rio Hondo (Oaxaca; 16° 10' 
N, 96°27' W), easternmost from Ocosingo (Chiapas; 
16°54' N, 92°06' W), and westernmost from Patzcuaro 
(Michoacan; 19°31'N, 101°36' W). 

Flying squirrels have been found at elevations 
from 840 m in Acahuizotla (Guerrero) to 2,800 m in 
Tepozan (Hidalgo) and 3,048 m in Cofre de Perote 
(Veracruz). They have a disjunct distribution, and 
populations are confined to the following eight, more 
or less well-defined, isolated mountain ranges: 1) Sierra 
Madre Oriental; II) San Luis Potosi; III) Sierra Gorda; 
IV) Eje Neovolcanico; V) Sierra Madre del Sur de 
Guerrero; VI) Sierra Madre del Sur de Oaxaca; VII) 
Sierra Madre de Chiapas; and VIII) Macizo Central de 
Chiapas (Fig. 1). Such mountain ranges vary in size, 
ranging from 801 to 67,142 km 2 , and are isolated from 
each other by belts of xeric (arid) or mesic (tropical) 
vegetation. The number of populations in each moun¬ 
tain range was weakly but positively correlated with 
their size (r = 0.87), and most localities were in the 
larger mountain ranges (VI - Sierra Madre de Oaxaca, 
and IV - Eje Neovolcanico). 

The distribution of flying squirrels is restricted to 
temperate plant communities, especially oak ( Quercus ) 
and mixed (oak-pine) forests (Table 2). All records 
were from oak (31%), pine-oak (56%), cloud (10%), 
and mixed (pine-fir-oak; 3%) forests. These forests 
have an annual precipitation between 600 and 1,200 
mm and average annual temperatures ranging from 12 


to 20°C. In Pinal de Amoles (Queretaro), flying squir¬ 
rels were found exclusively in oak forests and were 
recorded in 11 species of trees (five oaks and six pines; 
Table 2). The most abundant tree was the “escobillo” 
oak ( Quercus laurina), which accounted for 55% of all 
individuals sampled. A comparison between areas with 
and without flying squirrels showed that they differed 
in plant structure and plant species composition (Table 
3). The diversity of the species of trees was signifi¬ 
cantly higher in the areas without squirrels (t = -2.563, 
n = 11 ,p = 0.02). Habitats possessing flying squirrels 
were dominated by 0, laurina (84% of all trees); in 
contrast, in areas without squirrels, the dominant spe¬ 
cies were “roble” oaks ( Quercus crassifolia, ; 26.4%), 
“lacio” pines ( Pinus patula , 24%), and “prieto” pines 
( Pinus montezumae, 20%). Additionally, tree height, 
DPA, density and diversity of herbaceous plants, trees 
with cavities, and total number of cavities had lower 
values in areas without squirrels (Table 3). There was 
a positive significant classification of both sites with 
and without squirrels (X 2 = 27.89, d.f. = 4, P < 0.0000), 
based on DPA tree values, and greater density and 
diversity of herbaceous vegetation. 

Few flying squirrels nests have been recorded 
in Mexico, and all have been recorded in cavities of 
either oak or pine trees (see also Ceballos and Galindo 
1983; Ceballos and Miranda 1985). At Chapa de Mota, 
Mexico, all nests (n = 20) were found in dead pine trees 
in a pine-oak forest. In Pinal de Amoles, there were 
more trees with cavities and more cavities in plots 
with flying squirrels (Fig. 2, Table 3). The highest 
percentage of hollows was found in escobillo oaks (O. 
laurina , 74%), followed by roble oaks (O. crassifolia , 
16%), and black oaks ( Q. rugosa, 1.2%). Most (92%) 
nests were located in escobillo oaks and the remaining 
in roble oaks. All nests were lined exclusively with 
moss and possessed an average of two individuals. At 
Pinal de Amoles, the presence of oak trees seems to be 
a limiting factor for the distribution of G. volans , as we 
found nests only in the hollows of a single species of 
oak. The forest of escobillo oak presented the largest 
availability of hollows, as hollows naturally are formed 
in this species of oak. The number of hollows, in rela¬ 
tion with the total number of trees, is very small in the 


Table 1. Localities where flying squirrels (Glaucomys volans^) have been recorded in Mexico. Localities are organized by mountain range (Regions I-VIII) 
and indicate the primary plant communities and elevation. For forest type, the dominant tree group is listedfirst. Abbreviations are as follows: Latitude 
(Lat), longitude (Long), and elevation (Elev). 
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Region V (Sierra Madre del Sur de Guerrero) 

Guerrero 18) Acahuizotla 17°23' 99°27' Pine-oak 900 Goodwin 1961 

19) Agua de Obispo 17°17'99°22' Pine-oak, oak 960 Goodwin 1961 

20) Coapongo 17°23' 99°34' Pine-oak, oak 2,200 Goodwin 1961 
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Table 2. Habitat characteristics at 16 sites (100 m X50 m each) near Pinal de Amoles, Queretaro where Glaucomys 
volans were recorded as either present or absent. Flying squirrels were found in sites differing in tree species com¬ 
position and abundance when compared with adjacent sites without flying squirrels. Percentages are provided in 
parentheses. 


Tree Species 

Glaucomys Present 

Glaucomys Absent 

Tree Abundance 

Alnus jorullensis 

6(1) 

1 (<1) 

9(1) 

Arbutus glandulosa 

1(<1) 

44 (5) 

51(2) 

Cercocarpus macrophyllus 

0 

9(1) 

9(1) 

Pinus montezumae 

14(2) 

13(2) 

37(1) 

Pinus nelsonii 

7(1) 

180 (20) 

204 (6) 

Pinus patula 

46 (7) 

217 (24) 

523 (16) 

Quercus crassifolia 

18(3) 

237 (27) 

455 (14) 

Quercus laurina 

572 (84) 

95(11) 

1,812(55) 

Quercus mexicana 

1 C< 1) 

54 (6) 

56 (2) 

Quercus praineana 

2 (< 1) 

1(<1) 

7(1) 

Quercus rugosa 

144 (3) 

11(2) 

42 (5) 

Total 

678 

893 

3,207 


Table 3. Habitat characteristics of 16 sites (100 mX50 m each) near Pinal de Amoles, Quere¬ 
taro where Glaucomys volans were recorded as either present or absent. Abbreviations are as 
follows: H’ (Shannon Weiner Diversity Index), Hmax (maximum diversity), and J’ (evenness 
of species abundance). 


Characteristics 

Glaucomys Present 

Glaucomys Absent 

Number of tree species 

10 

11 

Tree density 

169.5 

223.2 

Number of individual trees 

678 

893 

H’ (Hmax) 

0.3 (1) 

0.8(1.04) 

J’ 

0.303 

0.771 

Average tree height (m) 

15.9 

12.1 

Diameter at breast height (cm) 

35.9 

24.2 

Number of hollows 

50 

3 

Shrub density 

193 

116 

Number of individuals 

770 

464 

Herb density 

117 

182 

Number of individuals 

468 

726 

Slope 

32.2 

39.8 
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Figure 2. Number of cavities suitable for flying squirrel nests in various species of 
oak trees at Pinal de Amoles, Queretaro. The different bars indicate one (diagonal 
hatched to the right), two (diagonal hatched to the left), three (hatched horizontally) 
or four (open) cavities per tree. 


area of Pinal de Amoles and is probably due to differ¬ 
ent natural and anthropological factors. For example, 
sawmills prefer trees of 30 cm dbh and 25 m of height, 
thereby limiting the number of trees that reach maturity 
and large size. Consequently, the number of den sites 
are limited, which in turn directly limits the utilization 
of habitat by flying squirrels (Muul 1968, 1974). 

Genetic diversity .—DNA sequences were ob¬ 
tained from 23 samples of G. volans representing the 
four populations examined in central and southern 
Mexico (Appendix I). Due to poor quality and quantity 
of some DNA samples, amplified sequences ranged in 
size from 243 to 1,165 bp. As a consequence of inser¬ 
tion/deletions, that were not a product of incomplete 
sequences, gaps were inserted to align sequences and 
resulted in a total alignment of 1,181 bp (including a 
portion of the proline tRNA region). 


Genetic distance values for selected populations 
are listed in Table 4. Values for within population 
comparisons ranged from 0.18% (Floresta, Chiapas) 
to 2.8% (Chapa de Mota, Mexico), and interpopula¬ 
tion comparisons ranged from 3.0% (Pinal de Amoles, 
Queretaro and Chapa de Mota, Mexico) to 5.5% (Pinal 
de Amoles, Queretaro and Floresta, Chiapas). 

The statistical parsimony network (Fig. 3) indi¬ 
cated that the haplotype shared by seven individuals 
from Floresta, Chiapas appeared to be the most basal 
haplotype based on the “star-like” topology with other 
haplotypes radiating from it. Based on the number of 
substitutions from the Floresta, Chiapas haplotype, 
samples from Chapa de Mota, Mexico appear to be 
unique and more divergent than samples from Pinal de 
Amoles, Queretaro or San Bartolo Morelos, Mexico. 
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Table 4. Average genetic distances estimated using the Tamura-Nei model of evolution (Tamura and Nei 1993) for 
selected comparisons of populations of Glaucomys. Ranges are indicated in parentheses, where appropriate. 


Comparison 

Average Genetic Distance 

Within populations 

Floresta, Chiapas 

0.18% (0.00-1.11%) 

Pinal de Amoles, Queretaro 

0.42% (0.00-0.83%) 

Chapa de Mota, Mexico 

2.81% (1.40-4.20%) 

Between populations 

Floresta, Chiapas - Pinal de Amoles, Queretaro 

2.56% (0.75-8.19%) 

Floresta, Chiapas - Chapa de Mota, Mexico 

4.10% (1.02-8.41%) 

Floresta, Chiapas - San Bartolo Morelos, Mexico 

1.70% (1.36-2.34%) 

Pinal de Amoles, Queretaro - Chapa de Mota, Mexico 

2.92% (0.94-4.30%) 

Pinal de Amoles, Queretaro - San Bartolo Morelos, Mexico 

1.12% (0.74-1.52%) 

San Bartolo Morelos, Mexico - Chapa de Mota, Mexico 

1.86% (1.53-2.68%) 

Chapa de Mota, Mexico - Oklahoma 

5.80% (5.48-6.10%) 

Oklahoma - All samples from Mexico 

5.40% (4.10-7.30%) 

Oklahoma - Nova Scotia 

5.20% (NA) 

All samples from Mexico - Nova Scotia 

10.10% (NA) 


^JK78479 Chapa de Mota, Mexico^ 



<^TK78478 Chapa de Mota, Mexico^)-©—©—©—©—o 


<^TK78476 Chapa de Mota, Mexico)—e—e—e—e—e—e- 


C^TK78480 Chapa de Mota, Mexico) o o o o o o 
<^TK90538 San Bartolo Morelos, Mexico^) 


TK78477 Chapa de Mota, Mexico^} 

-e—©—©—o 


TK78462, TK78455, TK78458, 
TK78459, TK78451, TK78456, 
TK78460 Floresta, Chiapas 


-e—C^TK78454 Floresta, Chiapas^) 


-©—<^TK90530 Floresta, Chiapas)> 


TK78450 Floresta, Chiapas^) 


TK90534 Pinal de Amoles, Queretaro^) 



TK90536, TK90535 Pinal de 
Amoles, Queretaro; TK78461 
Floresta, Chiapas 



TK90539 Pinal de Amoles, QueretartT^ 


TK90541 Pinal de Amoles, QueretarcT^ 


Figure 3. Network constructed using statistical parsimony methods. Individuals contained within the 
rectangle represent the most basal grouping. Large ellipses depict groupings of identical haplotypes. 
Small circles depict the number of base changes between adjacent haplotypes. 
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Discussion 


Geographic distribution and habitat use. —The 
dispersion route obtained for G. volans suggests that 
they invaded Mexico and Central America from the 
eastern United States, following the Sierra Madre 
Oriental through the center and south of Mexico to 
Central America (Ceballos et al. 2010). In Mexico, the 
subspecies of G. volans are endangered because of their 
restricted distribution, as they are found only in oak 
and oak-pine forest along the mountain range (Cebal¬ 
los and Navarro 1991). We support this conclusion, 
as our results show that the species has disappeared in 
several regions as a result of anthropogenic activities, 
especially deforestation. The natural distribution of 
G. volans is limited primarily to availability of cavi¬ 
ties. There is a relationship between the availability 
of cavities and the density of squirrels (Sawyer and 
Rose 1985). This situation could explain the scarcity 
of flying squirrels in most localities in Mexico, where 
the deforestation and changes in land use have reduced 
the surface cover of forests, eliminating trees that can 
be used as a refuge for this species. 

Genetic Diversity. —Based on the statistical par¬ 
simony network (Fig. 3), it appears that the ancestral 
haplotype may have originated in southern Mexico 
(Chiapas). Samples from Chapa de Mota, Mexico, 
were the most divergent based on the number of mu¬ 
tational steps separating them from other haplotypes. 
This divergence suggests that this population is more 
fragmented and has experienced less gene flow than 
populations from Queretaro, Morelos, and Chiapas. 
The samples from Floresta (with the exception of 
TK78461) grouped together as do the samples from 
Pinal de Amoles, Queretaro. 

Examination of intrapopulation variation (Table 
4) revealed that the population from Floresta, Chiapas, 
possessed the least amount of genetic divergence (aver¬ 
age genetic divergence, AGD = 0.18%), whereas the 
population from Chapa de Mota, Mexico, was the most 
diverse (AGD = 2.8%). Interpopulational comparisons 
showed the Chapa de Mota, Mexico, population to be 
the most different from the Floresta, Chiapas, popula¬ 
tion (average genetic distance = 4.1%). Although this 
study is limited by small sample sizes and geographic 
representation, it appears that substantial genetic 
variation exists among populations in southern Mexico. 


Levels of genetic divergence approach that between 
samples in Canada and the central United States (AGD 
= 5.2%) and between the United States and Mexico 
(AGD = 5.4%). Intrapopulation levels ranged from 
0.18% (Floresta, Chiapas) to 2.8% (Chapa de Mota, 
Mexico), indicating either a historical bottleneck in 
Chiapas or more widespread gene flow across central 
Mexico. Both hypotheses have their merits, and it is 
difficult to favor one over the other. First, the Floresta, 
Chiapas, population is located south of the Isthmus of 
Tehuantepec, long considered to have been isolated 
from southern Mexico (Ceballos et al. 2010). Second, 
the high levels of genetic divergence in the Pinal de 
Amoles, Queretaro, population may have been influ¬ 
enced by gene flow across the Sierra Madre Oriental 
and Transverse Volcanic Zone. No evidence of isola¬ 
tion by distance was detected. 

These findings support those of Kerhoulas and 
Arbogast (2010) who depict a separate evolutionary 
trajectory for North American and Mesoamerican 
populations of Glaucomys. Kerhoulas and Arbogast 
(2010) proposed that populations of Glaucomys from 
Chiapas, Mexico, are genetically divergent and perhaps 
isolated from other populations in Mexico, especially 
those from the Sierra Madre del Sur mountain ranges. 
Our data suggest that in addition to being divergent 
from populations in southwestern Mexico, populations 
from Chiapas are genetically different from populations 
located in the Trans-Mexican Volcanic Belt (Table 4). 

Conservation Status. —Currently, there are four 
protected zones with flying squirrel populations: the 
ecological reserve of Huitepec, in San Cristobal de las 
Casas; El Sepulcro, southeast to Ocozocuautla, both in 
the state of Chiapas; the Benito Juarez National Park in 
Oaxaca; and the ecological state park of Omiltemi, in 
Chilpancingo, Guerrero. However, these reserves rep¬ 
resent only a small portion of the geographic distribu¬ 
tion of G. volans. The current presence of flying squir¬ 
rels has not been confirmed in other reserves such as 
El Cielo biosphere reserve in Tamaulipas, where large 
tracks of cloud forests occur and flying squirrels are 
documented from subfossil remains (Martin 1955). 

A feasible strategy to maintain viable popula¬ 
tions of flying squirrels in Mexico will be to designate 
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protected areas throughout its geographic range. This 
especially is urgent in Chiapas, where suitable habitat 
is being destroyed at a rapid rate. Increasing the size 
of the Huitepec, Chiapas, reserve or other oak forests in 
the state is fundamental to maintain current populations. 
Important areas for the conservation of the species, be¬ 
cause of the extent of the remaining habitat, are located 


in central Mexico, in the states of Mexico, Michoacan, 
and Hidalgo. The selection of new reserves for flying 
squirrels can be combined with initiatives already es¬ 
tablished to complement the national protected areas 
system in Mexico (Ceballos et al. 1998; Chavez and 
Ceballos 1998; Ceballos 2007). 
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Appendix I 

Specimens examined .—Some samples used in this study were obtained from specimens that subsequently 
were released, and therefore voucher and museum catalogue numbers are not available for all specimens. All 
tissue samples were deposited in the Natural Science Research Laboratory, Museum of Texas Tech University, 
and are referenced with a museum identification number (TK). All specimens are from Mexico unless otherwise 
indicated. GenBank accession numbers are provided in parentheses following the TK number. 

Locality 6 in Table 1 and Figure 1 (N = 5).—Queretaro; Pinal de Amoles (TK90534, FJ376438; TK90535, 
FJ376440; TK90536, FJ376439; TK90539, FJ376437; TK90541, FJ376441). 

Locality 12 in Table 1 and Figure 1 (N = 5).—Mexico; Chapa de Mota (TK78476, FJ376456; TK78477, 
FJ376454; TK78478, FJ376455; TK78479, FJ376457; TK78480, FJ376458). 

Locality 13 in Table 1 and Figure 1 (N = 1).—Mexico; San Bartolo Morelos (TK90538, FJ376453). 

Locality 36 in Table 1 and Figure 1 (N = 11).—Chiapas; Municipio Comitan, 2 km S Floresta 16°28'32" 
N, 92°18'0" W, 2140 m (TK78450, FJ376444; TK78451, FJ376450; TK78454, FJ376445; TK78455, FJ376447; 
TK78456, FJ376451; TK78458, FJ376448; TK78459, FJ376449; TK78460, FJ376452; TK78461, FJ376446; 
TK78462, FJ376442; TK90530, FJ376443). 

Reference localities and outgroup samples .—UNITED STATES: Oklahoma; Seminole (TK26069, FJ376459); 
CANADA: Nova Scotia (AY928943; AY928973). 
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